Can disorder induce a finite thermal conductivity in ID lattices? 
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We study heat conduction in one dimensional mass disordered harmonic and anharmonic lattices. 
It is found that the thermal conductivity k of the disordered anharmonic lattice is finite at low 
temperature, whereas it diverges as k ~ TV ' 43 at high temperature. Moreover, we demonstrate that 
a unique nonequilibrium stationary state in the disordered harmonic lattice does not exist at all. 



Can disorder induce a finite thermal conductivity in 
one dimensional (ID) lattices? This question arose im- 
mediately after Anderson's finding of localization |l| . It 
has been commonly believed that disorder scatter normal 
modes and induce a diffusive energy transport that leads 
to a normal heat conduction. However, the early numer- 
ical as well as theoretical studies shown that the thermal 
conductivity in ID (mass) disordered harmonic lattice is 
proportional to A 1 / 2 , where N is the length of the lattice 
. The reason for this divergent thermal conductivity 
is not clear up to now. In fact, the existing theories ||]4|] 
do not prove that a unique nonequilibrium stationary 
state could be reached in a disordered harmonic lattice. 
Intended to obtain a finite thermal conductivity, Payton 
et al Q tried to add an anharmonicity to the disordered 
lattice. It was found that the anharmonicity did enhance 
the heat current, but Payton et al. were not able to ver- 
ify the validity of the Fourier law in such systems due 
to limited computer facility. Therefore, whether disorder 
can induce a finite thermal conductivity is still a puzzle. 



In past three decades, many works have been devoted 
to seek such a model whose heat conduction behavior 
obeys the Fourier law, i.e. J = —ttdT/dx. The primary 
motivation is to find out whether it is possible to prove 
and/or disprove a given ID many-body Hamiltonian sys- 
tem having or not having a finite thermal conductivity. 
If it is possible, then how can one calculate transport 
coefficients from microscopic Hamiltonian. The emer- 
gence of chaos in past years sheds some lights on this 
old but rather fundamental problem. It stimulated fur- 
ther study of heat conduction in non-integrable systems, 
aims at finding dynamical origin of the thermodynam- 
ical properties, such as irreversibility etc. ||. For this 
purpose, a wide range of non-integrable systems such 
as the Lorenz gas model H, the ding-a-ling and alike 
models 0, the Fermi-Pasta-Ulam (FPU) model g, the 
Frenkel-Kontoroval model H, the diatomic Toda lattice 
p0| , the Heisenberg spin lattice p| , the 4 lattice [12] 
and the sinh-Gordon and bounded single- well models [IS] 
etc, have been invoked. Most recently, the study has 



been extended to systems having periodic inter-particlc 
potential 0]. However, as we have already pointed out 
]l2| that non-integrability is only a necessary condition 
for the formation of a temperature gradient but not a 
sufficient one for a finite thermal conductivity. 

A rigorous proof of the necessary and sufficient con- 
dition for the Fourier law is still lacking, even though 
some progress in this direction has been achieved [ jl5| . 
Facing up such a situation, we have to rely on large 
scale molecular dynamics simulations to obtain a deep 
insight of underlying mechanism of the divergent and/or 
finite thermal conductivity. Recently p^ j, in searching 
for the underlying mechanism of the energy transport, 
we found that the (mass) uniform lattice systems can 
be classified into three categories: integrable one, non- 
integrable one with on-site potential, and non-integrable 
one without on-site potential. We discovered that the on- 
site potential plays a very important role. It destroys the 
momentum conservation and produces a diffusive energy 
transport, and thus leads to a finite thermal conductiv- 
ity. In non-integrable lattices without on-site potential, 
the scattering of solitary waves is found to be responsible 
for the divergent thermal conductivity k ~ N 0A3 in the 
FPU and alike models. 

In this Letter, we would like to study the effect of disor- 
der on heat conduction in ID lattices. The Hamiltonian 
of the models we are considering can be written as 



H, 



(i) 



where V(xi-\,Xi) stands for the interaction potential of 
the nearest-neighbor particles. In our study the mass of 
particles is given by: 



mo + A(i?, - 0.5), 



(2) 



where A is a parameter adjusting the amplitude of the 
random mass, and R t is a random number distributed 
uniformly in the interval of (0, 1). 

Disordered anharmonic lattices. In this case the inter- 
particle potential takes the form as that in the FPU 



f 



model, 



1 



V(xi-i,Xi) = -(x 



2 + ?(* 



(3) 



/3 = 1 throughout this paper except otherwise stated. 
We call the model disordered FPU (DFPU) model. In 
our numerical simulations the Nose-Hoover thermostats 
[Il6| were put on the first and the last particles, keep- 
ing them at temperature T+ and T_, respectively. The 
eighth-order Runge-Kutta algorithm was used to solve 
the coupled differential equations. All computations were 
carried out in double precision. Usually the stationary 
state set in after 10 6 ~ 10 7 time units, thereafter the 
time average (Jj) (Jj = iidV / dxi + i) was found to be 
site independent, and thus is denoted by J. 
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FIG. 1. JN versus N (a) and temperature profiles (b) for 
DFPU model (A = 0.4) and the FPU model at high temper- 
ature, T+ = 3 and T_ = 2. The solid lines in (a) are the 
best fitting ones, both have same slope a — 0.43 ± 0.01. (c) 
JN versus N for the DFPU model with different values of 
disorder at low temperature, T+ = 10~ 3 and T_ = 5 x 10 -4 . 
The lines are drawn to guide eye. (d) Same as (b) but for the 
case of low temperature. 

In Fig. |l|(a), we present JN versus N for the FPU 
and the DFPU models with T+ = 3 and T_ = 2. Both 
models have a similar (almost indistinguishable) temper- 
ature profile (Fig. 0(b)). In both cases, the temperature 
gradient is proportional to -/V -1 , thus JN oc k. Best fit- 
ting in FigJ|(a) yields k ~ Ar°- 43±0 01 . This implies that 
in both models the thermal conductivity diverges in the 
same way. More numerics show that in this high tem- 
perature regime, the disorder does not help converge the 
thermal conductivity, it only decreases the value of heat 
current. 

Interesting things turn out when the temperature of 
the thermostats is decreased. In Fig.^J(c), JN versus N 
is plotted for the DFPU model at three different val- 
ues of disorder A = 0.2, 0.3, and 0.4, respectively, with 



TV 



10 _d and T_ 



5 x 10" 



When N < N c , 



JN saturates, namely k becomes an ^-independent con- 
stant. This result indicates that the heat conduction 
obeys Fourier law when the size of the ID lattice becomes 
larger than N c . However, in this temperature regime, the 
FPU model and the harmonic lattice have the same heat 
conduction behavior . This is reflected in the temper- 
ature profile shown in Fig. 0(d). It demonstrates that no 
temperature gradient can be formed, and the stationary 
state corresponds to T = (T+ + T_)/2. 




JN increases monotonically with N. When N > N c , 



FIG. 2. J(N)/J(2N) (for N = 600) versus the average 
temperature T = (T+ + T_ )/2 in the semi-logarithmic scale 
for the DFPU (solid circle) (A = 0.4) and the FPU (solid tri- 
angle) models. The FPU model shows a transition to a har- 
monic lattice with the decrease of the temperature. However, 
the DFPU model changes from a divergent thermal conductor 
to a normal one as the temperature is decreased. 

In FigJ|, we show J(N)/J(2N) versus T = (T + +T_)/2 
for the FPU (solid triangle) and the DPFU (solid circle) 
models. Where J(N) is the heat current for the lattice 
of length N. The transition of the FPU model to a har- 
monic lattice happens when the temperature is decreased 
below a threshold value of T c ( w 10" 2 ). This is illus- 
trated in Fig. 2 by an abrupt decrease of J(N)/ J(2N) 
from 1.5 to 1. On the other hand, within almost the same 
temperature range the heat conduction in the DFPU 
model undergoes another transition. This transition is 
demonstrated by an abrupt increase of J(N) / J(2N) from 
about 1.5 to 2 when T < T c . This implies that the DFPU 
model transforms from an abnormal thermal conductor 
to a normal one. 

The above results of the FPU model can be understood 
from the scattering mechanism of solitary waves. As we 
discovered |l2j that the scattering between solitary waves 
from opposite directions causes an energy loss that gives 
rise to a temperature gradient. However, the momen- 
tum conservation in the FPU model prohibits a diffusive 
energy transport. The interaction of the solitary waves 
depends on the initial excitation. When T < T c , the 
solitary waves cannot be excited any more, the energy 



2 



is transported by linear excitations - phonons, thus the 
FPU model behaves like a harmonic one. 




FIG. 3. Energy distribution in space at time t = 200 
for high excitation - solitary waves (a) and low excitation - 
phonons (b) for a DFPU lattice with disorder A = 0.4. (c) 
The energy of solitary waves on different site. The dotted 
lines are the best fitting ones with the data in initial stage 
and in final stage, respectively. The different slopes (given in 
the figure) show that the energy decays in different rate for 
solitary waves (initial stage) and phonons (final st 



In order to clarify the underlying mechanism of the 
heat conduction in the disordered model, we study the 
dynamics of a single excitation. Figure ||(a) shows a snap- 
shot of energy density Hi in a DFPU lattice at t = 200 for 
an initial excitation on the left end, i.e. pi — 3,Pi = 0, 
for i 7^ 1. A solitary wave is excited and propagates to 
% = 250. If we decrease the energy of initial excitation be- 
low a threshold value, we find that the energy is no longer 
localized, namely the solitary wave disappears. Instead, 
a traveling wave packet, linear excitations, is created and 
moves diffusively (see Fig. ||(b)). However, unlike in the 
case of the FPU model, when a single solitary wave or a 
linear wave packet travels along the disordered lattice, its 
energy will decrease due to the scattering from the dis- 
order. An interesting phenomenon is that solitary waves 
and phonons obey different scattering laws, as is shown 
in Fig. ||(c). In this figure we plot the energy versus the 
distance of the solitary wave excited in Fig||(a). We see 
that both in initial and final stages, the energy decays 
in power law H s ~ N~ a but with different exponents, 
a ss 0.15 in the initial stage and a « 1 in the final stage. 
The former represents the scattering of solitary waves 
by disorder, and the latter represents the scattering of 
phonons by disorder as in the final stage the energy be- 
comes so low that the solitary waves cannot be excited 
any more. 

With above picture, we can explain the numerical re- 
sults of the DFPU model shown in Figs. [l]-|| At high 
temperature (Fig. |l](a-b) and Fig. ||), energy is trans- 
ported mainly by solitary waves which encounter two 



kinds of scattering: scattering from solitary waves and 
scattering from disorder. The scaling law of energy loss 
caused by the two scattering are N~ a but with differ- 
ent values of a, a ~ 0.5 for the former |l2), a « 0.15 
for the latter (see Fig 3(c)). In the thermodynamic limit 
the main cause of energy loss comes from the scattering 
of solitary waves. Thus at high temperature the DFPU 
model shows the same scaling law as the FPU model as 
is shown in Fig |l](a). In this case the disorder cannot in- 
duce a finite thermal conductivity, it decreases the total 
heat current only. 

The situation becomes different at low temperature 
(Fig. [l] (c-d) and Fig.||). In this case, the solitary wave 
is hardly excited. The energy is transported by phonons. 
The scattering of phonons by disorder becomes a dom- 
inant factor for energy loss. Thus the energy decays as 
iV" 1 as is shown in Fig. ||(c), which implies a finite ther- 
mal conductivity. This prediction is confirmed by a direct 
calculation of JN given in Fig. 0(c). It is natural to re- 
late the saturation phenomenon of JN in Fig. |l|(c) to 
the localization of phonons in the system, wherein the 
localization length I is an important scale. If the length 
of the lattice N is shorter than I the thermal conductiv- 
ity is expected to be divergent. On the contrary, if the 
lattice is much longer than I then a normal heat con- 
duction could be realized. This is indeed the case shown 
in Fig.|l|(c). The threshold value N c is approximately 
500, 700, and 1000 for the three different values of the 
disorder A = 0.4,0.3, and 0.2, respectively. Extensive 
numerical simulations confirm that N c oc 1/A. This is of 
order I, because I oc (m)/y/ ((m — {m)) 2 ) ~ 1/A |Q. 

Disordered harmonic lattices. It was shown theoret- 
ically ppj that the thermal conductivity diverges as 
k ~ N ^ in ID disordered harmonic lattice (/3 = in 
Eq.([j])). The proof was based on an assumption that a 
unique nonequilibrium stationary state could be reached 
in such systems. However, to the best of our knowledge, 
this assumption has not been verified numerically. This 
motivated us to do further investigation. 

We find that a stationary state, with erratic temper- 
ature fluctuation, can be set up as is shown in Fig. 
||(a). The states starting from the same initial condition 
(T° = 0.01) are almost indistinguishable at different time 
t = 10 6 and t = 10 7 . (Please note the inset for magnifi- 
cation.) This guarantees the existence of the stationary 
state. Unfortunately, such a state depends sensitively on 
initial condition . In the same figure, we show another 
stationary state (dashed line) formed from a different ini- 
tial condition T° = 3. It differs greatly from that one of 
T° = 0.01. This result means that the nonequilibrium 
stationary state in the disordered harmonic lattice is not 
unique. 

To obtain a clear picture, we turn to a simple model, 
a lattice consists of only two types of atoms, one with 
mass m = 1 the other with mass m = 0.5. A seg- 
ment of light atoms is embedded in the middle of other 
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two segments of heavy atoms. Two ends of the lattice 
are put in the thermostats with the same temperature 
T + = T- = 4. Two stationary states from two differ- 
ent initial conditions T° = 0.1 and T° = 1 are shown 
in Fig.|](b). Again, no unique stationary state can be 
formed. However, adding a small fraction of anharmonic- 
ity, f3 = 10~ 8 , to this lattice produces a complete differ- 
ent picture. With this tiny anharmonicity, we can obtain 
not only a smooth temperature profile but also a unique 
stationary state, see Fig. ||(b). 
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FIG. 4. (a) Temperature distribution of a disordered har- 
monic lattice (A = 1) from different initial conditions, the 
solid lines for T° — 0.01 at 10 6 and 10 7 time units (see also the 
inset for the magnification of a small interval), respectively, 
and the dotted line for T° = 3 at 10 7 time units, (b) The 
temperature profile (solid line) of a harmonic lattice consists 
of three uniform segments. The left and the right segments 
are atoms of the mass m = 1 while the middle one consists of 
atoms of mass m — 0.5. T+ = T_ = 4. The dashed line is for 
the same system but with a tiny anharmonicity /3 = 10 -8 . 

The above results show that although the anharmonic- 
ity alone is not enough to yield a finite thermal con- 
ductivity, it plays a crucial role in establishing a unique 
nonequilibrium stationary state in disordered lattices. 

In summary, we have studied heat conduction in ID 
disordered harmonic and anharmonic lattices. Our nu- 
merical results show that at low temperature disorder can 
induce a finite thermal conductivity. The magnitude of 
the disorder does not affect the results in the thermody- 



namical limit. It only determines the localization length 
and threshold length of a lattice having a finite thermal 
conductivity. However, at high temperature, the disor- 
dered anharmonic lattice shows a divergent thermal con- 
ductivity, which is similar to that of the FPU model. In 
addition, we provide a numerical evidence showing that 
the nonequilibrium stationary state in a disordered har- 
monic lattice is not unique. 
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